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ABSTRACT 


Both integrator and canceller types of sampled 
analog comb filters implenented by a canonical second 
order recursive er SU Er have been studied 
experimentally and theoretically. Two 96 stage CTD 
(charge transfer device) delay lines are used in the 
experimental study. Theoretically, an analysis based 
on a modification of the digital recursive filter 
theory was developed and accounted for the major 


features of the experinental results. 


Three comb filter applications have been 
demonstrated. The first two used the integrator type 
comb filter for the improvement of the signal-to-noise 
ratio of a periodic signal contaminated by noise and 
for the sorting of a periodic signal of selected PRF 
from other periodic signals. The third application 
used a canceller type comb filter to eliminate a 


Sinusoidal interference from a desirable signal. 
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DEFINITION OF TERNS 


Several equations are presented throughout this 
thesis. Many of the terms used in these equations are 
generally represented in most texts and papers as Greek 
symbols. To simplify the typing requirements, these terms 
have been redefined. The following list shows the standard 
Symbols which are used throughout this text. 


PI = 3.141592654 


f = frequency = 1/T 
w = 2*PI*f 

e = ezp 

5 = jw 

4 = ехр(5Т) 


6 
MHz = Mega Hertz = 10 Hertz 


3 
K = Kilo Hertz = 10 Hertz 


db = decibels 

Pf = Pico Farrads = 10 ^ Farrads 
t = Tao 

7 z Any Arbitrary Starting Tine 
CTD = Charge Transfer Device 


SNR = Signal To Noise ratio 
CIT ےک‎ 5112 0-2-2350 


Transfer Function = V IN. 
output input 
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A. GENERAL 


One of the major considerations in svsten design is 
signal processing. For all tvpes of svstems, it is 
necessary to consider how the signal will be represented 
(analog, digital, or Sampled analog), received and 
conditioned to minimize errors and noise; and finally, how 
the desired signal will be recognized and/for discriminated. 
All these factors must be considered whether the system 
under concern 1S a communication, radar, sonar, accoustical, 
Optical, or some other type. Although the signal processing 
techniques vary greatly in all the previously mentioned 
systems, the operations performed have a great deal in 
common. The many processing operations can be grouped into 


four catagories: 


Two types of filters: 


1, Nonrecursive 
a. Convolvers 


b. Correlators 
2. Recursive 
Two types of spectral analysis: 
1. Pourier type of transforas 


2. Other tvpes of transforms 
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This thesis is concerned with the sampled analog 


recursive filter. Specifically, tkə recursive comb filter 
is studied both experimentally and theoretically. The 
feasibility of this filter's application is also 
demonstrated. 


HE RECURSIVE COMB FILTERS 


A comb filter is a class of filters whose transfer 
function is periodic in the frequency domain. There are two 
classifications of comb filters: canceller or bandstop, and 
the integrator or bandpass. The canceller is characterized 
as having a very narrow but a high attenuation region 
between neighboring teeth and good transmission within the 
passband region. The integrator is the opposite of this in 
that the passband region exhibits good signal passage in 
peciodically spaced narrow bands of frequencies and a 
relatively wide attenuation band in between neighboring 
teeth. Comb filters have been developed for both analog and 
digital filtering systems. The development hd these 
circuits has been based on the feedforward, frequency 
sampling, fast Fourier transform, and recursion technique of 
implementation. The common feature of 120111. 
implementaticn for all these devices is that they all 


require delav lines as the basis of their operation. 


The delay devices which are used vary with the systens 
encountered. For example, shift Tegisters are used in 
digital systems whereas sampled analoq systems use the 
bucket brigade device (B8D), charge coupled device (CCD), or 
the serial analog delay (SAD) device. Due to the sampling 
technique of these sampled analog devices, they exhibit both 
sampled analog properties and some of the digital properties 


as well. 
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The comb filter was analyzed for a sampled analog filter 
to determine its frequency response characteristics. The 
technique used in designing this type of filter vas the 
Bilinear Z transform which is discussed in Chapter 2. This 
filter was implemented as a recursive filter which is a 
filter that is characterized as a device whose output 
depends on the present input, past output, and its present 
output. It can also be ¿implemented as a non-recursive 
feedforward filter and used for spectral analysis through 


the chirp-Z-transforn (CZT). 


E PURPOSE OF THE STUDY 


During the past two years, recursive comb filters have 


been studied at th2 Naval Postgraduate school. Vi 
[10] .. . ۱ 

Tamsa-ad first studiəd the use of charge transfer device 

delay lines for the implementation of sampled analoq 


DT 
recursive filters. | B. Fruend studied the first order 
recursive filter and emphasized the comb feature. ورا‎ 

(8) ۱ 
Saetre extended the study to the second order recursive 


comb NeT and demonstrated the application of the 
6] 

canceller type comb filter to MTI radar. S. Holmes 

developed a computer aided design procedure for the 

recursive comb filters and studied their sensitivity to 

device errors and variations. The status of tha recursive 

comb filter study before the beginning of this thesis can be 

Summarized as follows: ~ 

-Theory- 
The digital recursive filter theory was modified to analyze 


the sampled analog recursive comb filters. While the 


(ә 





modified theory accounted for the comb feature of the 
filter, it was not able to explain three major features of 
sampled analog recursive comb filters which are different 


from an ideal comb filter. 


xThe peak of the frequency characteristics do not remain 
constant as in an ideal comb filter. They become frequency 


dependent in the sampled analog implementation. 


*The frequency characteristics do not remain zero between 
neighboring teeth as in the ideal case. Deterioration of 


these nulls generally increases with frequency. 


*The Shape of the comb teeth do not remain constant as in 


the ideal comb filter. These also change with frequency. 
-experimental study and Application- 


Only the canceller type comb filter was experimentally 
studied in detail and its application to MTI radar was 
demonstrated. The Integrator comb filter was not studied 
thoroughly. 

Therefore, there are three objectives for this thesis: 


1. Improve the theoretical analysis. 


2. Experimentally study the integrator comb filter and 


demonstrate its applications. 


3.Demonstrate other applications of the canceller comb 


filter and discuss some of its applications. 


D. OUTLINE 


Chapter 1 presents the background and purpose of this 


thesis. Chapter 2 summarizes the basic ideas about digital 
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filters together with a discussion of the mapping transforms 
(H(s) into H(z)) vhich are used for the design of recursive 
conb filters. In Chapter 3, the design procedures for both 
the canceller type and the integrator type recursive comb 
filters are discussed. The experimental results with comb 
filters are presented and compared with the theoretical 
analysis in Chapter 4. Chapters 5, 6, and 7 present 
applications o£ comb filters. Chapter 5 specifically 
demonstrates the use of the integrator comb filter asa 
means to improve signal-to-noise ratio of a pulse train 
Signal which has been contaminated with white Gaussian 
noise. The integrator comb filter is used as a pulse train 
sorter to select a particular pulse train from among varied 
frequency pulse trains. A third application of the comb 
filter is presented for a canceller comb filter in Chapter 7 
where the device is used to eliminate sinusoidal 
interference from a second, desired sinusoid. Finally, 
Chapter 8 summarizes the results and offers suggestions for 


further studies. 
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A. INTRODUCTION 


The theory of digital recursive filters has been vell 
developed in several text books and published papərs as vell 
as some past theses written here at the Naval Postgraduate 
School(i.e. see referencəs 6,7,and 8). Por this reason the 
theory of recursive filters will not be covered in detail, 
although a general overviev of sampled analog theory vill be 
presented beginning with the transform theory. This will 
lead into a discussion of the transforms wnich are used for 


designing both digital and sampled analog filters. 


B. BASIC TRANSFORM THEORY 


The analysis of analog circuits is highly dependent on 
the Laplace Transform to ease the computations which are 
involved as well as to simplify the circuit synthesis of the 
desired equations. For these same reasons transforms have 
been developed for use in digital circuits. It has been 
found that these same transforms can be used in sampled 
analog circuits. The Z Transforms vhich have been developed 
can be considered as steppingstones from either the Laplace 
S-plane domain, the discrete time domain  (t-nT), or the 
continuous time domain into the Z-plane domain (see figure 
Zol). 
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Sampling 


5:6 


rs Bilinear 
Stamdard 


Figure 2.1 Transform Steppingstones 


ESR TRANS FORMS 


There are a number of transforms which can be used to 
change a transfer function in the S-domain to a transfer 
function in the Z donain: 

1. Impulse Invariant Transformation (Standard Z) 

2. Bilinear Z Transformation 

3. Mapping of Differentials 

u. Matched Z Transform 

5. Others | 
The most popular of these are the Impulse Invariant and the 
Bilinear Z transforms shown in figure 2.1. These mapping 


functions are given below: 





The Standard Z Transform: 
ን L) 
Ehe Bilinear 2 Transforn: 


2) 1 - 2-1 ( 
۲" ) 1 + 2-1 ነ 


D: TRANSFORM PROPERTIES 


The inherent properties of these types of transforms are 
similar although they do have some pronounced differences 
which make one technique more desirable than the other for 
some applications. These properties of both transforns are 


shovn in Tables 1 and 2. 


This thesis utilizes the Bilinear Z Transform technique 
because it maps the entire left half S-plane into the entire 
unit circle and there is no aliasing effect. These were the 
two major reasons for using this transform, other reasons 


are shown in Tables 1 and 2. 
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TABLE 1 


TRANSFORM እጅ .”።፦:። || PROPERTIES O | PROPERTIES 


Bilinear Z Simple algebraic mapping between 


continuous and digital filters. 


Maps entire jw axis in the S-plane to 


the unit circle in the Z plane. 


Realizable, stable continuous circuits 
map into realizable, stabie digital 


systems. 


Videband sharp cutoff continuous filter 
can be mapped into videband sharp 
cutoff digital filter without 


allasing. 


There is a non-linear relation between 
the analog and digital frequencies 


resulting in frequency warping. 


The amplitude response of the analog 


systems being transformed must be 
piece-wise constant, otherwise the 
digital frequency response will be 
warped. (See ref. 3, pg. 222-223) 


Neither the impulse nor the phase 
response of the analog filter is 


retained in a digital filter. 
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TABLE 2 


TRANSFORM سس‎ መመ . PROPERTIES 


Standard Z Encounters an inherent aliasing problem 


unless the sampling frequency is high. 


The input signal nust be band linited 
betveen -PI/T Hz and PI/T dz for the 
frequency responses of the analog 
filter and its digital filter design 


equivalent to correspond. 


The impulse response of the digital 
filter is a sampled version of its 


analog counterpart. 


To use this transfornation, tne order 
of the numerator must be less than the 
order of the denominator of the analog 
filter to avoid intolerable aliasing 


problems in the digital filte 
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-1 
E. PHYSICAL SIGNIFICANCE OF Z 


-1 
The operation of Z is an important property in 


designing a sampled data circuit. Mathematically, Z is 


Simply a complex variable; but, physically, it represents a 
unit of delay for a sampled data circuit. How then can such 
a function be implemented? There are two techniques which 
are used to implement this delay physically. A digital 
circuit uses a shift register or memory. For sampled analog 
Circuits the family of Charge Transfer Devices (CTD) delay 


lines are used. 


F. DIGITAL FILTERS VERSUS SAMPLED-ANALOG FILTERS 


There are major differences between digital filters and 
sampled-analog filters. However, much of the digital filter 
design theory is carried over into sampled analog filters. 
A digital filter receives an analog signal, quantizes it 
(i.e. represents the signal amplitude bv a series of pulses 
which make up a coded binary number corresponding to the 
analog amplitude) and then perforus its filtering action. 
This quantizing action is accomplished by an 
analog-to-digital (A-D) converter. It is generally an 
operation which introduces errors and extraneous electrical 
noise into the filter system. A sampled-analog filter 
eliminates this A-D conversion by sampling the signal's 
amplitude and storing a quantity of charge proportional to 
the signal amplitude. The storage (or delay action) is 


accomplished by the charge transfer devices for analog 
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filters and by shift registers for digital filters. 


The Trelationship betveen digital filters and sampled- 
analog filters that are interchangeable is seən in the 
design techniques. The well developed digital filter design 
theory is applicable to sampled- analog design as well with 


some modification. 


Ga cos UNIT CIRCLE 


When mapping from the S-plane to the Z-plane, it is 
important to consider which regions in the S-plane and 
Z-plane will correspond. This knowledge is important to 
consider the stab: litv of the filter designed according to a 


transfer function in the Z domain. 


For the Bilinear Z transforn, the entire left half 
S-plane maps into the unit circle in the Z-plane, shown in 
figure 2.2, whereas the right half S-plane maps outside the 
unit circle. 





S-Plane Z-Plane 


FIGURE 2.2 Bilinear Z Mapping 
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This is an inportant result, since the mapping of the 
entire jw axis onto the perineter of a unit circle preserves 
the frequency selective properties of the analoq systen. 
The mapping of the left half S-plane inside the unit circle 
in the Z-plane ensures that stable analog filters will nap 


into stable discrete time filters. 


The desirable “mapping characteristics described above 
are not inherent with the Standard Z transform, shown in 
figure 2.3. The Standard Z transform only maps periodic 
strips of the left half S-plane into the unit circle in the 
Z plane. The riçht half of the same strip maps outside of 
the unit circle. This is the reason for the possibility of 
aliasing at improper sampling frequencies. These strips are 
2PI/T Hz in width and the first desirable mapping strip is 
between -jPI/T and jPI/T radians in the S-plane. 


iW 


Imaginary 





S-Plane Z-Plane 


Figure 2.3 Standard Z Mapping 
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III. FILTER DESIGN 


A. SAMPLED ANALOG FILTER DESIGN 


i? 1 2 0 0116 1 8 


The design and implementation of a filter can be 
approached fron many aspects. Basicallv, the numerous 


techniques which are used fall into one of two catagories; 


the Direct Method or the Indirect Method. Rabiner and 
[3] 


Gold in their excellent textbook have discussed these 
methods in detail. The procedures used in this thesis were 
a version of the indirect method using the previously 
discussed transforms. Figure 3.1 shovs the procedural flov 
of these design techniques. In this section, the design and 
characteristics of the recursive comb filter 1 be 


covered. 
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DIRECT TECHNIQUE [INDIRECT TECHNIQUE 


SPECIFICATION IN THE 
H(s) FREQUENCY DOMAIN 


| 
| 
| 
FREQ, | 
一 | APPROXIMATION IN THE 
| S-DOMAIN His) 
| USING THE FOLLONING 
| FILTER PROTOTYPES: 
| *BUTTERNORTH 
| *CHEBYSHEV 
| 


şehir 
BESSER 


MODEL THE ABOVE 
ANALOG FREQUENCY 
DOMAIN FILTER CHARACTERISTIC 
DIRECTLY INTO THE 
Z-DOMAIN 


SAMPLED ANALOG DIGITAL FILTER ACTIVE FILTERS 


FILTER H(s) H(z) VARIOUS CIRCUIT 
[7] STANDARD-Z CONFIGURATION USING 
BILINEAR- OP, AMPS, 
VERY SIMILAR E B 
10 DIGITAL FILTERS 
ETC, 
IMPLEMENTING | CANONICAL REALIZATION: 
HARDWARE: PARALLEL, CASCADE 
RETICON 
BBD 


Fr 
Ku, 





Figure 3.1 Design Methods 
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A filter vhich depends on its present input, past 
output, and also its present output is defined as a 
Recursive Filter. This definition implies that there are at 
least tvo operations vhich nust be performed to synthesize 
such a device; feedback and delay. In addition to these two 
operations, a Recursive Filter also includes multiplication 
and addition. All these operations are evident by studying 


the Z transfer function as seen in equation 1 below: 


32. trade 
Hz unu 070 07” 1 
9 1200 D. 2 > 


This equation is developed fron either the analog 
lov pass or the analog high pass transfer function using the 
Bilinear Z transformation technique. The derivation of this 
equation is a common one and is found in most text books. 
However, for completeness the reader is encouraged to read 


Appendix B of reference 7 pages 99 through 120. Freund, B. 
7 

e J in his thesis has derived, using the Bilinear Z 

transform method, the transfer function for the first and 

second order lov and high pass filter. Since the derivation 

of equation 1 has been thoroughly explained by Fruend the 

design techniques here were used to determine the 


coefficiənts for the proper response that vas desired. 


The relationship betveen the nunerator coefficients 
is an important one since the numerator polynomial 
determines the location of the zeros of the transfer 


response. In the derivation of H(z) the normalized transfer 
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function shoved that the relationship betveen Se a , and a 


9551. 2, 1, respectively. Às shown below, this gives a 
double zero on the unit circle either at Z214jO or 221-0 


depending upon the sign of the 2. coefficient. 


For جج‎ 


ዘ 
ፎጋ 


-1 
1:22 +Z 


2 
Multipiy by Z 


2 
2 *22*1 = 0 
OLr 


2 +1) (5 +1) = 0 
27” 


therefore 


sSiönyilarıly fər ከ 
2 
Z -2Z+1 = 0 
(2-1) (2-1) = 0 


4 = 1, Z = 1 
1 2 


Plotting these values of Z on the unit circle reveals that 


there is a dcuble zero at Z= +1 or at Z = -1. This 


indicates that either a low or high pass filter is 
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determined simply by the siqn of 27 (plus for lov pass and 
minus for high pass). See figure 3.2. Other nunerical 
relationships betveen ... a and E will result in the 


zeros becoming complex, and/or shifting to either inside or 


outside the unit circle resulting in deterioration of the 


transfer function stop bands. 


Z-Plane Z-Plane 
=] 1 =| | 
/ 
a. Tov Pass Filter Do 1121 255 er 
aq = +0 a, “ -4.0 
2,= 2-4 il-Za = + 4 


FIGURE 3.2 Zero location for low or high pass filter. 


The "b" coefficients also effect the transfer 
Characteristics, although much more differently than the a" 
coefficients. Тһе denominator coefficients affect the shape 
Of the pass band region width. By either increasing or 
decreasing the pass band, the 3db cutoff frequency is also 
changed. This change in the pass band shape is accomplished 
by the shifting of the pole locations within the unit circle 


as seen by the calculation below: 
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Denominator equatioa 


-1 -2 
1۲95 2 +b Z - 0 
1 2 


+ 
Multiplying throuqh by 5 
2 
2 *Ы 24Ы =0 
1 2 


Factoring 


b, + jVü4b, - bi bb - jV4b, - b 


( Z вт) (а, ہچ رکچ‎ 2 


2 2 2 


Thus the problem of the filter design becomes a 
problem of selecting the proper value of the denominator 


coefficients in the Z transfer function. 


Insight in determining what the coefficients should 
be for a particular response can be obtained from the 
comparison of the unit circle pole-zero plot and the desired 
frequency response of the filter. This is done by first 
placing the zeros on the unit circle at the appropriate 
frequency. Usually the zeros are desired on the unit circle 
to afford as nearly as perfect a stop band as possible. The 
poles are placed inside the unit circle at the proper 
frequency and their distance from the center determines the 
sharpness (quality factor) of the peaks (or pass bands). 
Figure 3.3a below shows a desired low pass frequency 


response, However, when some insight concerning the pole and 
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zero location is desired, the procedure is as follovs: 


1. Place the zeros at tne proper position on tne 
unit circle for either the high or lov pass 


filter charactəristic. 


2. Locate poles inside the unit circle at the 
frequency at vhich they occur on the frequency 


response. 


Imaginarv 





ЈЈЈЛ 


a Desired Response bD- Estimating the Pole ama 
Zero Location” 


FIGURE 3.3 Graphically Estimating Pole Zero Locations 
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This procedure can be performed ir reverse also. 
For example,if the poles-zero locations are known, the 
frequency response can be estimated by moving 
counter-clockvise around the circle starting at Z= ۰ 
The amplitude is increased as a pole is approached and 
decreased as a zero is approached. Midway between the poles 
and zeros, an arbitrary median amplitude value is 
Maintained. If the pole is closer to the perimeter of the 
circle the peaks are Sharp (high Q) and if they are closer 
to the center of the circle the peaks are not sharp but 


wider (lov 0). 


4. Integrator and Canceller Recursive Comb Filters 


— q == — ہے ےو‎ ይ. መጻ m e Q= sama sm سم سن = .€—.—,—.—..-...- انس وص‎ ss HELD ge eee =መ፦ 


The Recursive Comb filter can be configured as 
either an integrator or a canceller type of filter. The 
difference between these types of filters is simply in the 
relationship between their stop band and pass band widths. 
The integrator is characterized as having a very narrow pass 
band (thus implying a very high Q circuit) anda very wide 
stop band. The canceller is the opposite of this with aà 
very vide pass band and a narrow stop band. These 


characteristic responses are shown in figure 3.4. 


5 








Canceller 





۲۶۴٢ 7 


Hiyunr.: Keocuürsive Comb Filter Types 
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5. Coefficient Determination 


One portion of this thesis deals with Integrator 
Recursive Comb Filters, and the demonstration of tneir use 
as a device to improve signal-to-noise ratio of a pulse 
train signal added to white Gaussian noise. The Z transfer 
function is given by equation 1 (given in Chapter 2) and the 
numerator coefficients are known (1, 2, 1 respectively). 
The design of the filter as discussed earlier is a problen 


of determining the denominator coefficients ә апд .. This 


is not a trivial problem since there are very many possible 


combinations of ር and 25 each between +1 and +2. Where 
then does one begin searching for the proper coefficients 


that will give an integrator Comb filter response and also 
have a very narrov pass band? Due to the massive amount of 
calculations required to plot a frequency response of 
equation 1 the problem lended itself well to to a computer 
iteration. The fortran program shown in Appendix B is a 


programing of equation 1. 


To obtain the fraquency response from equation 1, 
Z=Cos ( wTJ +jSin (wT) was substituted into the transfer 
function resulting in H(jw). This vas then programed and 
the only inputs that were required were the numerator anda 
denominator coefficients. The output is a plot of the 
frequency response on a linear graph and the amplitude is 
normalized to the maximum amplitude encountered throughout 
the frequencv range of interest. This program will 
automatically plot the rəsults on the textronix terminal of 


the CP/CMS time sharing systen. 


As noted earlier equation 1 is the transfer 
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enarseteristic of a digital filter design technique. This 
equation applies equally well to the Sampled-Analog filter 
design technique, although a slight modification must be 
made to account fər the sample and hold circuitry which is 
used in the filter to reproduce the filtered signal at the 
output. In. the frequency domain the sample and ከ014 


characteristic is a SinX/X function which affects the filter 


response in a multiplicative manner. This characteristic 
has been included in the fortran progran. One other 
characteristic which was necessary 1600) model the 


Sampled-analog filter was the adjustment of the frequency of 
the sin and cos arguements to give the proper number of comb 
teeth, corresponding to the number of delays which were used 
in each stage of the filter. A second order filter was 
designed with 96 delays per stage. This resulted in a comb 
response with 96 teeth between D.C. and the sampling 
frequency. The 96 teeth are a result of the repetitive 
nature of the bilinear design technique. This 
characteristic vas also modeled into the fortran program 


which is shown in Appendix B. 
6. Filter Design 


The computer program of Appendix B will give the 
frequency response of the desired Z transfer function, but 
as previously mentioned the coefficients are required input 
data. Thus, a method is needed to determine vhere to start 


searching for the proper coefficients. This problem was 


6) 


( 
answered by the result of a thesis by Stacy 7۰ Holmes 


In his research he derived a means to analize Sampled-Analog 


filters based on the coefficients d and SC In his 
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dissertation he introduced a b-plane triangle of stability 


(6) 


(shown here in figure 3.5 ın. From this graph a rough 


estimation of the denominator coefficients can be obtained. 


Although the coefficients which are obtained in this manner 
are very rough figures it does provide a starting point to 
use the previously discussed computer model. Several "b" 
coefficients can be used in the program until the desired 
frequency response is obtained, prior to building the 
filter. Since this study is of recursive Comb Integrators 


and Cancellers, from the graph the = coefficients are 
limited to values of approximatly 40.2 to -2.0 where as the 
> coefficient can very fron 41.0 to -1.0. This technique 
is ,needless to say, time consuming and tedious; however, 


the b-plane triangle of stability decreases the coefficient 


search immensely. 
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EN 72 ۱ 
_ > 
Dividing Isocontour between LS 
Cancellers (above) and 
Integrators (belov). 





Figure 5 5053 b-Plane Triangle of Stability 
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5317 COMPUTER AIDED FILTER DESIGN 


One of the intentions in this study is to design a 
Recursive Integrator Comb filter with a very narrow passband 
and a broad stop band and demonstrate its use as a device 
for improving Signal-to-Noise ratio of a pulse train input 
Signal (a topic of a later chapter). Suffice it to Say that 
the type of filter design desired can be found with 


coefficients that are located in the region below the .05£ 
5 


isocontour on the b-plane triangle of stability shown in 


6 
figure 3.5 | Several values of coefficients vere tried in 


the computer program. these coefficients are as shovn in the 
table belov: 


TABLE 3 
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Figures 3.6 through 3.25 show the computer results for 
the coefficients listei above. These plots give an 
indication of the results which can be expected 
experimentally from a filter with these coefficients and 


this number of delays. 
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Figure 3.6 Comb Frequency Spectrum (Theoretical) 
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Figure 3.9 Comb Frequency Spectrum (Theoretical) 
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Figure 3.10 Comb Frequency Spectrum (Theoretical) 


45 





Normalized Amplitude 


Оут 


0 


= 
123 eo 3e 
Frequency (KHz) 


a Akbar 


Figure 3.i1 Comb Frequency Spectrum (Theoretical) 
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Figure 3.12 Comb Frequency Spectrum (Theoretical) 
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Figure 3.14 Comb Frequency Spectrum (Theoretical) 
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Figure 3.15 Comb Frequency Spectrum (Theoretical) 
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Figure 3.16 Comb Frequencv Spectrum (Theoretical) 
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3.18 Comb Frequency Spectrum (Theoretical) 
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Figure 3.19 Comb Frequency Spectrum(Theoretical) 
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IV. EXPERIMENTAL FILTER STUDY 


À. CIRCUIT CONFIGURATION 


The Recursive Comb filter vas connected as shown in 
Figure 4.1 below. 


Summer Summer 


Input 





Figure КЕТ ет Circuit Commicirat ion 


wer circuit 6522 10112 1 15 the conical ሮ ርን ተ 
realizationş using the Bilinear Z transform of either the 
high or low pass S-Plane analog filter transfer function. 
The sample and hold circuit was necesserv to give better 
nulls in the frequency response characteristics of the 


filter. If this circuit had not been used, the sumation of 
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the signals at the operational amplifier summer vould be the 
addition of a strict analog signal and a sampled analog 
signal. This introduces errors and rapidly deteriorates the 
nulls of the frequency response. With the sample and hold 
circuit, tvo similar signals are added in the operational 
amplifier. To further reduce errors, the clock frequency of 
the sample and hold circuit was adjusted to be the same as 
the sampling frequency of the delay elenents. The 
coefficients vere composed of carbon resistors with five 
percent tolerences. The resistors were arranged as voltage 
divider networks where the values were adjusted bv combining 
resistors in series and measuring the tap voltage (which was 
a portion of the entire voltage aeross tue series resistance 
combinations). The tap voltage was determined by 
multiplving the total voltage on the entire series 
combinations by the desired coefficient value, and after 
this value vas measured at one of the tops, the entire 
combination vas placed into the circuit. It was necessary 
to obtain the coefficients in this manner, instead of using 
variable potentiometers for coefficients, since these latter 


devices were found to be frequency dependent (as determined 


577 
1 ) 


in earlier studies 


The delay device was the Reticon SAD-100 serial analog 


delay. Jn the circuit diagram of Figure 4.1, the total 


n 
delay is represented as Z . Since the Reticon SAD-100 has 


ninety -six bits of delay, the value of n in this circuit is 


ninety-six clock periods for each delay stage. 
BZ EXPERIMENTAL DATA 
The circuit which has been described above was studied 
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to determine its frequency response characteristics and the 
dependence of these responses on the denominator 
coefficients. The numerator coefficients were Ot no 
interest here since it was determined in Chapter Two tnat 
their relationship must be fixed to give the comb response; 


and, that changing the sign of 22 defines either a high or 
low pass filter (+ for low pass and - for high pass filter) 


characteristic. Therefore, for all data presented here, a 


low pass filter was used. 


The data plan is the same as the plan used for the 
computer design section (table 3). This vas done to compare 
the results of the computer program. Table 3 has been 


repeated here for convenience. 


TABLE 3 





The experimental data is shown on the following pages in 


Figures 4.2 through 4.17. The two sets of curves for b, 2 


መሙ 
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+1 were not taken since the device vas unstable with these 


coefficients. Accompanying each frequency plot is a plot of 
the poles and zeros of the transfer function for the 
coefficients used in the respective response. This was 
included to obtain a better understanding of the reason for 
the shape of the filter teeth. The Z plane plot starts at 
D.C. at Z= 1 + jo and one counterclockwise revolution, 
around the circle, increases the frequency to fs/96 , or one 
complete cycle of the frequency response that it 
accompanies. The clock frequency for the following figures 
vas 400 K Hz. 
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Figure 4.2 Comb Frequency Spectrum (Experimental) 


and Pole Location 
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Figure 4.4 Comb Frequency Spectrum (Experimental) 
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Figure 4.5 Comb Frequency Spectrum (Experimental) 


and Pole Location 
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Figure 4.6 Comb Frequency Spectrun (Experimental) 


and Pole Location 
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Figure 4.7 Comb Frequency Spectrum (Experimental) 


and Pole Location 
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Figure 4.8 Comb Frequencv Spectrum (Experimental) 
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Figure 4.11 Comb Frequency Spectrum (Experimental) 


and Pole Location 
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Figure 4.12 Comb Frequency Spectrum (Experimental) 


and Pole Location 
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Figure 4.13 Comb Frequency Spectrum (Experimental) 


and Pole Location 
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Figure 4.14 Comb Frequency Spectrum (Experimental) 


and Pole Location 
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Figure 4.15 Comb Frequency Spectrum (Experimental) 


and Pole Location 
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Figure 4.16 Comb Frequency Spectrum (Experimental) 


and Pole Location 
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Figure 4.17 Comb Frequency Spectrum (Experimental) 


and Pole Location 
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1. ata Analyses 


It appears , from the data in the previous figures, 


that changing the sign of > causes the fılter response to 
Change from a canceller to an integrator type fiiter. This 


is not completely true. The reason that the filter is 
changing back and forth between canceller and integrator can 
be explained by referring to Figure 3.5 of Chapter 3. In 
this figure, there is an isocontour denumerated .25fs which 
distinguishes between the two types of filters under concern 


here. The fact that ” is held numerically close to this 
isocontour causes the filter characteristic to jump back and 
forth across this isocontour. However, when Ba = 0.2 3 e 


is greater than 0.6 in absolute value, the filter response 
remains as an integrator since these coefficients fall below 
(integrator side) the .25fs isocontour. There is a large 
difference in the shape of the frequencv response teeth just 


dascribed. The changing of the sign of b: from minas to 
plus causes the poles to change from being real(minus Ra 
coefficients) to being imaginary (plus b coefficients) even 
though both curves are integrator responses. 

Another point which was noted was, that as the value 
o£ P decreases negatively, the slopes of the curve skirts 
become  steeper. This was expected, since, as the b-plane 
triangle of stabilitv indicates, when b becomes more 


- 


negative, (holding e constant) the value of the isocontour 
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ratio becomes smaller. This indicates that the three 
decibel cutoff frequency is decreasing (f was held constant 
s 


throughout the experiment). This, in effect, causes the 
pass band to decrease, bringing about a sharper response 
curve. This can be seen by comparing Figures 4.3, 4.4, 


4.6, and 4.8 and their counterparts which are plotted with 


b lus. 
2 p 


Figures 4.14 and 4.16 show the nulls beginning to 
deteriorate quite rapidly. This is due to the shifting of 
ehe Zeros around the unit circle from the Z = +1 + jo point. 
Since there is a real pole relativelv close to the perimeter 
Oft the unit circle in these two figures, the pole 
immediately becomes effective; allowing sone of the signals 
at these frequencies to pass through. However, as the 
frequency was increased still further, the second zero is 
encountered and a second null obtained. The deterioration 
became increasingly apparent as the frequency was increased. 
An example of this deterioration is shown in ۳ 1 1 ۱۲ ۵ ۵۵ 
This figure shows the entire spectrum of the filter response 
from DC to fs = 400 K hertz. Also shown in Fiqure 4.18 is 
the deterioration of the nulls due to the shifting of the 


zeros and some inefficiency of the Reticon filter. 


The shifting of the zeros may have been due to stray 
capacitances caused by the connecting leads of the circuit. 
This would cause the coefficients of both the transter 
function numerator and denominator to vary and, in the case 
of the numerator coefficients, cause the zero locations to 
shift off the reai axis. This results in an imaginary term 


in the value of Z, causing a shift around the unit circle. 
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277 PILTER SELECTİON 


As will be discussed in the following chapter, an 
integrator recursive comb filter is desired vith very narrov 
pass band teeth to demonstrate the Signal-to-Noise ratio 
improvement ability of the comb filter. The type of 
response chosen was the filter with E = -0.2 and Ss = “0:6. 


This filter has a very narrow three decibel pass band 
(approximately 417 hertz vide) and a vide stop band 
(approximately 3.91 K hertz wide). This is shown in Figure 
4.6. Figure 4.19 shows the entire frequency response from 
DC to 400 K Hertz. The zeros are relatively unaffected at 
the higher frequencies although the amplitudes of the teeth 
are increasing between 80 K hertz and 110 K hertz. This vas 
possibly due to the frequency dependence of the b 
coefficients since, as explained earlier, these coefficients 
affect the teeth shape. However, for the purposes of the 
next chapter, this was acceptable. The discrepency in the 
teeth amplitude at approxinately 200 K hertz is due to 


instrumentation errors only. 
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V. APPLICATIONS OF IHE INTEGRATOR TYPE RECURSIVE COMB 
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A. INTRODUCTION 


Signal-to-Noise ratio (SNR) improvement is a major 
concern in radar and communication design problems. In this 
section an attempt has been made to use a Charge Transfer 
Device (CTD) to improve the SNR. The input signal that was 
chosen was a pulse train. A pulse train was chosen since it 
is commonly used in radar applications as well as the fact 
that it has a discrete frequency spectrum with many 
frequency components. This latter fact demonstrates the 
intent of this section more explicitly, as will be seen 


later. 
B. FOURIER ANALYSIS AND FILTER CONSIDERATIONS 


The result of a Fourier analysis of a pulse train is 
given by equation (1) below and its development is shown in 


Appendix A. 


Vt 1... t t/2) د1‎ Ыб6/2 
k 1 A kw. t/ 2 


where: 
t =pulse width 
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I= pulse repition frequency (PRP) 
V= pulse amplitude 
w= 2PI/T 


The line Spectrum is shown in Figure 5.1. Note that the 
amplitude is a function of the familiar SinX/X curve (where 
X = kwt/2 for this explanation ) and that all the even and 
odd harmonics are present since the pulse train in this 
situation is neither an aven nor an odd function of tine. 
Further note that the value of ር is directly related to the 


pulse width (t) and inversely proportional to the pulse 
period T. I£ the pulse period (T) increases, the band width 
2 : pecterun ( band width is from 0 hertz to  1/T) does 
not change, but the number of spectral lines between freq.- 
O and freq.- 1/T increases. If this pulse train vere 
transmitted through a very noisy environment it would be 
ሀ. ےہ‎ Cult, ٩ not impossible, to detect it unless 
Special filtering techniques were employed. However, if the 
frequency spectrum of the pulse train could be matched with 
the frequency spectrum of a filter which will attenuate all 
other frequency components in between the desired signal 
components, then the system noise will have minimum effect 
on the transmitted/received signal. The problem then is to 
find a filter which has the same frequency spectrum 
۰۰۶ teristics as that of Figure 5.1. A filter with this 
type o£ spectrum is the Integrator Recursive Cənb filter. 
This filter is periodic in the frequencv domain and is 
capable of nearly perfect attenuation vithin the stop band. 
Since the spectrum of the pulse train nas discrete frequency 
Components over a very vide frequency range a very narrov 
comb filter pass band response is desired to capture the 
signal's spectral components and at the same time attenuate 
as much of the undesired spectral components as possible. 
The characteristic response of the Integrator Comb filter is 


such a response. The characteristic response of the 
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Integrator Comb filter is shown in Figure 5.2. This 


response can be produced very easily using sampled analog 


techniques. Sampled analog techniques are highly desirable 
since there is no analog-to-digital or digital-to-analog 
conversion required. Also, the frequency spectrum of the 


comb filter can be shifted without changing the relative 
amplitude of the pass band, simply by changing the sampling 


frequency £ . Since the comb filter frequency response 5 
5 
a finite band width per tooth, some of the noise spectral 


components will pass through the device with the signal!s 
spectral components. The major portion of the noise should 
be attenuated since the pass band widths are much narrower 


than the stop band widths of the filter response. 
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Figure 5.2 Integrator Comb Filter Response 


CE DEVICE SELECTED 


The device selected to compose the filter was the 
Reticon SAD-100. This device jis a monolithic integrated 
circuit which provides storage of analog information. The 
device has 96 bits of delay and provides a first-in, 
first-out storage for these 96 samples. The specification 
sheet on the SAD-100 indicates that there are 100 bits of 
delay, two of which are used for input and output junctions 
leaving 98 storage cells. Earlier work with this device had 
indicated that there were 96 storage cells. The data which 
was obtained throughout this study proved this to be 


accurate. 


Some of the important characteristics of the SAD-100 are 


given below: 
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Maximum Delay at ü KHz sampling rate 25 msec 


© e o e € 


1577:7ü7 Samplıing Rate..................... 12 MHz 
Dynamic Range...... 25555. s . ONG 
NEM Clock Frequəncy...............1... 2.0 2 
"55:76 Cell Capacıtance..............a... he 2 pf 
HüDUc Capacıtance.......... 15: m 505 5700 pf 
5521060 Capacı tance- ........ -............. 75:75: pf 


The coefficients which vere used are: 


=1.Q‏ و 
0 
02ح ظط 2.0 = a‏ 
1 1 
a z1.0 9 ን.‏ 
2 2 
These coefficients give the most desirable frequency‏ 


response as shown in Figure 5.3. The stop bands are nearly 
perfect zeros and the pass bands very narrov vith very sharp 


slopes on both sides of their respective center frequencies. 


The filter vas operated at a sampling frecuency of 390 K 
hertz during the measurements. Thus, if the pulse train is 
to pass through the filter unattenuated and without 
distortion, the spectral components of the pulse train must 


fall directly on the pass bands of the filter response. 
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The pass bands of the Comb filter are periodic with 


respect .to the sampling frequency (f ) divided by the 
5 
number of delays in the SAD - 100. The filter"s first tooth 


center frequency is given by equation 2 belov: 


> = f /Number of Delays 2 
5 


This calculation gives a center frequency of the first 
tooth of the filter at . =4063 hertz. Since thə follovinq 


teeth are repetitive with respect to this frequency, the 
following teeth center frequencies are given by M*f , where 
i 1 


M is an integer equal to 1, 2, 3, ... 


For the signal to pass through the filter, the spectral 
components of the signal and those of the filter must 
coincide. For the pulse train to coencide with the filter 
frequency spectrum, it must have a pulse repitition 
frequency of 4063 Hertz or some integer multiple of this 
frequency. The frequency selected was 4063 Hertz and its 
frequency spectrum is shown in Figure 5.4, Figure 5.5, and 
Figure 5.6. Figure 5.5 and 5.6 are an expanded scale of the 
input signal , whereas figure 5.6 is the entire spectrum 
from DC to 450 KAZ. A comparison of figures 5.3 and 5.4 
will reveal that the spectrums are matched perfectly. The 
time domain representation of both the input and output 
pulses are shown in Fişure 5.7. In this figure only one 
input pulse is shown, and above it, one output pulse. The 
slight distortion in the output pulse is due to the loss of 
some of the higher frequencies resulting from an imperfect 
frequency match. This loss is caused by slight drifting of 


the pulse generator from the selected frequency, as well as 
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the Zeros of the stop band at the higher frequencies rising 
above zero. Figure 5.8 shovs a series of input and output 
pulses. Figure 5.9 demonstrates the signal attenuation 
when the input signal frequency spectrum is not matched to 
the filter spectrum. The input frequency in this case is a 
3.56 K hertz pulse train ( shown on the bottom of figure 


5.9); note that the output signal is completly Lost (shown 
in the top of figure 5.9). 


94 





اا یں UN ተ‏ کت دئاع | 

lA. SE T. ۳-9 ۱ ነ 
TT E eg muz ا‎ eri” 
Bex a a cup S ETT ту rro اص‎ 
A A ACA E 








əpn3TTduy 032 ۲ ۲ 6111 ۷ 


95 


44.09 48.70 


30.57 40. 5 


32.5 
Frequencv (KHz) 


28.44 


24.38 


.13 12.19 16.25 20.5 


8 


4.00 


0.0 


4063 Hz) 


irequencv Spectrun (PRF- 


Figure 5.4 Pulse Train l 





- e 
PE Ди 2” 
یہ‎ 


"ووو 


aw 





ከ ኤሌ emi rarn Frequency 


Spectrum (Experimental) 


227«8ሃ 





1043/10xHZ SÖkHZ 


EHiçür6 570 input Pulse Train Frequency 


Spectrum (Experimental) 


96 








T —>"-.F—S —_ 一 一 一 一 - — —‏ م ت ee —— — >” = > TEST‏ سس سم 


hai put Pulse 





purer- S Filter Input and Output 


97 








5171777” Gancellatizon of the Input Signal 
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The input noise which was used was produced by General 
Radio Company and was a type 1390-5 random noise generator. 
This device produces a white Gaussian noise over 3 frequency 
range of up to 5 M hertz. The frequency spectrum of the 
Gaussian noise was plotted on the X - Y recorder and is 
shown in Figure 5.10. This graph is plotted from Zero hertz 
to 500K hertz. The slight dip from approximatly 100K hertz 


to 300K hertz is due to the noise generator characteristics. 


The noise as seen in the frequency response of Figure 
5.10 was passed through the Reticon integrator Comb filter. 
As to be expected, the filter frequency response should be 
Sambar to tha -lof figure 5.3. The actual noise output 
frequency response of the filter is shown in Figure 5.11. 
This plot is from zero hertz to 50K hertz. A comparison of 
this and figure 2's response reveals that the base band has 
widened and that the peaks have become jagged. This is due 
to the noise having a very low input amplitude of 0.6 volts 
RMS. Thus, the product of the very low input Signal and the 
filter pass band tooth on a normalized basis causes the 
small amplitude skirt portion of the pass band to remain 
approximatly the same as it vas before but, the product of 
the very large center portion of the filter pass band with a 
very small input spectral components causes the large center 
frequency amptitude to decrease considerably giving the 
effect of widening the pass band of the filter on a 
normalized basis. Hovevər, this in effect is not the case, 
(Note that all of the actual curves plotted have been 


normalized to ease the comparisons which have been made.) 
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Figure 5,12 graphically illustrates the attenuation of 
the noise signal with the input noise spectrum and output 
noise Spectrum plotted on the same graph. The graph is 
plotted from zero to 50K hertz. The noise input spectrum is 
the top curve and the output spectrum is below it. This 
same input and output noise are shown in Figure 5.13 in the 
time domain. Note that although much of the noise is passed 
through the filter, much of it has been cancelled out. In 
Figure 5.13, the bottom is the input noise and the top, the 
output noise. 


/ 
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5:15 515 lime Domain of input and Output Noise 
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D. SIGNAL-TO-NOISE RATIO MEASUREMENTS 


It has been determined at this point that the Integrator 
Comb filter will act well in cancelling much of White 
Gaussian noise. But, can a signal be ietected well enough 
to determine its presence ; and, if so, what type of signal- 
to -noise ratio improvment can be obtained? These are very 


important questions and vill be pursued next. 
The frequency spectrun of the input pulse train and 


nəise vas added vith a resistive netvork as shovn in Figure 
5.14 below. 


1000 ohms 


Signal 


1000 ohms 
f م‎ Signal plus Noise 
Noise Output 


1000 ohms 


Figure 5.14 Network for adding Signal and Noise 
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The resistances were necessary to prevent the noise from 
feeding into the signal source as well as to prevent the 
signal from feeding into the noise source. The 1.0K ohm 
resistor to ground was necessary to prevent the relatively 
low input inpedance of the filter from changing too 
drastically. 


P700£c 0.15, Fıguüre 5.16, Figure 5.17,.and Figure 5.18 
are photographs of the signal plus noise input and output. 
The input Signal was decreased from 2 volts peak to peak 
(figure 5.15) to 0.1 volts peak to peak ( figure 5.18). 
This was done to give an indication of how well the filter 
will reproduce the signal even though it is buried in a 
relatively high noise environment. As seen in figure 5.18 
the output pulses ( bottom trace ) are much more noticable 


than the input pulses. 


Figure 5.19 demonstrates the cancellation of the input 
signal, which is added to the Gaussian noise, when the input 
frequency is not matched to that of the filter. Note that 
the output trace (bottom ) is only composed of filtered 
noise whereas the top trace ( input ) is composed of both 


the signal and noise. 


The previous data gives a visual indication that the 
signal is reproduced and also gives some indication of how 
well the signal is reproduced. However, a more precise 
method in answering the question of how well the signal is 
reproduced is to compare the signal's input and output 
frequency spectrum's. The following photographs show the 


frequency spectrum comparisons of the input and output. 
Figure 5.20 is the frequency spectrum plot of the 


output pulse train alone. It covers a frequency range from 


zero to 450K hertz. Note that there is a large impulse 
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390K hertz. This is due to the sampling clock PRP feeding 
through the filter; not a characteristic of the pulse train. 
Further note that the spectrum decreases to approximatly -67 
dB at its minimum value which occurs at 390.0K hertz. This 
is as predicted since the Fourier analysis (as seen in 
appendix A ) predicts that the spectrum will have a  SinX/X 
envelope with the first null occuring at the reciprocal of 
the pulse width ( i.e. 390.0K hertz ). Figure 5.21 shows an 
overlap of both the input and output spectrums. The two 
curves are very much alike with the exception that the 
output spectrum ( top curve ) is approximatly 10 dB higher 
than the input spectrum (bottom curve ). This is due to the 
filter having a gain of approximatly 10 dB. Figure 5.22 is 
the same output spectrum of figure 5.20 except that the plot 
extends only to approximately 20K hertz. Finally, Figure 
5.23 shows the input and output plots overlapping. Again 
note the increase in the peak amplitudes due to the gain of 
the filter. Thus this comparison of both the time domain 
and frequency domain proves that the pulse train is 


reproduced exactly. 
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Figure 5.16 Signal Plus Noise 
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Figure 5.18 Signal Plus Noise 
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Ergure 5.23 Output and Input Frequency Response Comparison 





E. FREQUENCY CONSIDERATION 


To this point there has been no comparison of signal 
plus noise in the frequency domain. The following data will 
show that the input signal is uneffected by a noisy 


environment and is reproduced exactly by the filter. 


Figure 5.24 15 a frequency spectrum of the input signal 
plus noise from zero hertz to 50 K hertz. The spectrun 
shown is a continious one since Gaussian noise has a 
continious spectrum; but , in this graph the pulse train 
spectral components are added to the nolse spectrum causing 
the spectral lines to protrude higher than the noise 
spectral amplitudes. For a signal - to - noise ratio 
improvement it vill be necessary for the Integrator Conb 
filter to attenuate the noise spectral components. This is 
in fact what has been done as seen in Figure 5.25. All the 
noise spectral components between the signal spectral lines 
have been attenuated to nearly zero. Note, however, that 
there is still noise passed through the filter since the 
pass hand teeth have some finite width. Also, the width of 
the output spectrum teeth appear to be much wider than the 
filter pass band teeth. This is due, as explained earlier, 
to the small signal levels. Thus from the previous data, 
the siqnal is in fact reproduced exactly since all the 
spectral components of the pulse train signal are passed 
through the filter but also along with this reproduced 
signal there is noise passed through; however, the noise has 


been reduced considerably. 
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P. SNR IMPROVEMENT USINS A SECOND ORDER FILTER 


It is evident at this point that there has been an 
improvement in signal - to - noise ratio, since the signal 
has been reproduced exactly and much of the noise has been 
attenuated. The question still remains as to how much 
improvement has been obtained. The answer to this question 
was obtained by making voltage measurements with an RMS 


voltmeter. The measurements were as follows: 


Pulse Train Voltage (RMS ) . 104 .880 
Noise Voltage ( RMS ) 2205 15 
Filter Noise Voltage < RMS ) = 1025 


Since these are voltage measurements the signal-to-noise 


ratio is calculated using equation 3 below. 


SNR(dB) = 205 (5378 V = LOG ¥ .- ) 3 
10 sig 10 noise 


Thus the values are as jiven below: 


SNR, = -16.12dB 
ın 


SNR 
out 


232 dB 


The difference of the output SNR and the input SNR is 


the inprovement. For these conditions it is: 


SNR = 13.79 dB 
improvement 


Although a 13.79 8۵8 SNR improvement is not a great feat 


the fact remains that, for a first generaticn application of 
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the Reticon SAD -~ 100 as an Integratcer Recursive Comb 
filter, it has performed well. Room for improvement in its 
operating techniques as well as the filter design techniques 


still remains. 


С. FIRST ORDER FILTER SNR IMPROVEMENT CURVES 


The work which has been done above vas all with a second 
order integrator filter. The following data was taken with 
the Reticon SAD-100 operating as a first order filter having 


numerator coefficients as: 1 “a... 1:0. Thə b 


coefficient was variable for this demonstration. The 


transfer function for this filter contiguration is shown 


below as equation 4. 


1 + 2-1 
H(z) شت‎ ü 
1 + 5,271 


The purpose of this section vas to obtain some insight 
of the SNR improvement ability of the first order filter. 
The experimental data is shown in Appendix C. This data vas 


plotted and is presented below in figure 5.26. 
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The graph of figure 5.26 above is plotted with the b 
coefficient on the abscisa as the independent variable and 


the SNR improvement as the dependent variable. the data was 


plotted in this manner since 62-19, coefficient determines 
the shape of the filter passband teeth as well as whether 


the filter is a canceller (t b's) or an integrator (- b's). 


Thus by varyinc Ss the filter will either pass more or less 
ncise component through with the pulse train spectral 


components giving a variation in the SNR improvement. 


Three curves are plotted each with varying SNR. 
input 


values. Although curves A and B only differ by 0.4 db, the 


input signal and noise voltages doubled in curve B giving 
approximatly the same input SNR. This vas done to see the 
effects the input amplitude had on the amount of 
improvement. Curve C is a plot with the input 8 
approximatly three times the value of either curves À or. 8. 
It vas interesting to note that all the curves are within 
2db or less of each other as seen in figure 5.27. This vas 
an expected result since the filter pass band testh are the 
same for the respective b coefficients Tregariləss of the 
input signal amplitude. The SNR improvement is much higher 
for negative b coefficients because of the filter teeth 
becoming much narrower (i.e. a higher Q integrator circuit) 
and attenuating more and more of the noise. Since the three 
curves of figure 5.26 are all relatively close to each 
other, it vas felt that the variations are due to 
instrumentation measurement errors and the fact that the 
noise is such a highly random signal that the measurement of 
the noise changes slightly each time a new measurement is 


made on the same signal, These small changes in 





measurements cause a noticable logarithmic variation due to 


the nature of the logarithmic conversion characteristic. 
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A. RADAR 


Now that the Integrator Recursive Comb filter has been 
demonstrated as a device with potential to improve 
Signal-to-Noise ratio, to what practical use could such a 
device be put? One application which has been discussed in 
the literature (Recent Advances in the Synthesis of Coapb 
filters) is their use for Radar application for video 
integrators. They are used in this application to 
accomplish linear addition (integration) of the radar 
signals which are consistant from sweep to sweep and at the 
same time the random signals ( noise ) will add in a root -- 
sun -square form. The net result is an improvement in the 


SHR as shown in the previous section. 


B. ELECTRONIC COUNTERMEASURES 


Another application which is proposed bv this report is 
the use of the Integrator Recursive Comb filter as a Pulse 
Sorter of Radar pulses tə speed up the acquisition of Radar 
Signals in an Electronic Countermeasures situation. This 
concept of a pulse Sərter has been denonstrated and the 
results are shown in Figure 6.1. The systems pulse trains 
added together make up the first four traces. There are four 


separate traces to show each pulse train as a stationary 
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signal on the scope. This was done by synchronizing the 
oscilloscope separately with each input signal and 
reexposing the same negative to the scope trace, after 
synchronizing and shifting the position of the stationary 
pulse train. After the first exposure, the filter output 
was moved off the scope szreen, the next pulse train vas 
synchronized to the szcope, (to keep it stationary) the 
entire trace was moved vertically down on the screen, and 
the film vas reexposed to the screen. Thus, the same 
negative was exposed four times. The frequencies of each 


pulse train from top to bottom are: M = 8.46K hertz, = = 
1000 hertz, z = 2000 hertz, and Se = 1500 hertz. The last 
trace is the filter output of the captured pulse train (5 
in this situation). Note the correspondence between the 
pulse locations in the top trace = and the output trace 
(bottom trace). This is the captured signal. The others do 


nət pass through the filter since they fall in the filter's 
stop bands and are attenuated. Thus, by using the Integrator 
Recursive Comb filter, a particular pulse was sorted fron 
among several and interference reduced from non-synchronous 


pulse train. 
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Erpure 60. imPulse Sorting Illustmation 
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Ideally, the Reticon filtering will not pass through 
the system at SEIL c But, there are some frequency 
relationships between the desired passed pulse train and the 
undesired signal so that the filter will allow for passage 
of some of the undesired components. The components that 
will pass are those of the signals that have a whole 
multiple spectral component of the filter's clock frequency. 


For example, the frequencies which were used above fu. z 
8.86 K hertz, > = 1000 hertz, =. z 2000 hertz, and H = 1.5 
erez) , £ , £ , and 7 have spectral components that are 


2 3 
د‎ = 8Khertz, ila = 8Khertz, and "i = 12Khertz and will 


fall within the pass band of the filter teeth. The effects 


of these components will cause some (as yet undetermined) 
distortion. The distortion was not noticable in the data 
taken in this thesis. An indepth study of this distortion 
would be most interesting as a thesis chapter in the future. 
Figure 6.2 below shows the spectral components that fall 
within the filter pass bands and very close to the captured 


signal Bh, 
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The previous figures have demonstrated that the 
filter can select a specific pulse train signal and pass it 
while attenuating the others. The signal vhich is selected 
by the filter is reproduced at the output; but, how 
accurately is the signal reproduced? Figure 6.3a shovs, on 
a linear plot, the line spectrun of the input pulse train. 
This is the signal that was eventually passed through the 
filter to give the output trace in Figure 6.1 above. 
Although the filter will pass multiples of signals starting 
at 406K /96 Hertz, the primary frequency selected for this 
demonstration was 8.46 K hertz. With this frequency, every 
other tooth of the comb filter frequency response will 
contain a spectral component of the 8.46 K Hertz input 
signal. This was done to ease the comparisons which are 
made later in this report in determining the filter's 
selectivity. Figure 6.3b shows the output of the filter 
using the input signal as shown in Figure 6. 3a. Note that 
there are some very small spectral components at 
approximately 4K Hertz, 12K Hertz, 20K Hertz, etc. These 
components are that portion of the comb filter's teeth pass 
bands that are not filled with spectral lines of the signal. 
Instead, there are sonə noise components (possibly filter 
qenerated) vhich fall into these unused teeth and are being 
passed through. In comparison of these two signals, Figure 
5.4 shovs an overlay of both the input and output line 
spectrums. It is important to note that the signals" 
spectral lines match perfectly and at the proper frequency. 
This shows that there is no frequency distortion as the 
signal passes through the filter, but does not give any 
indication of phase relationship. Note, however, that the 
spectral amplitudes are very different. This amplitude 


change is a result of the filtering system gain. 
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In summerizing the observations of Figures 6.1,6.2,6.3, and 
6.4 is noted that: first, the filter reproduces the input 
signal without frequencv distortion; second, the filtering 
system anplifies the signal, and third, the filtering system 
introduces some small amount of noise into the signal. NOW 
that it is known that the filter can reproduce the input 
pulse train with an added benefit of amplifying it; what 
will be the effects of two input pulse trains with equal 
PRF"s and amplitude be on the filter? The following data 
will answer this question. The spectral components of the 
8.46K hertz pulse train are shown in Figure 6.5a. It was 
expected that if a second pulse train of an equal PRF was 
injected at the filter input, the spectral components would 
simply add. This effect was demonstrated and is shown in 


Figure 6.5b. 





۳ 1۱۲ zu peecr ar OUtput With Two Input Signals 


129 





ای 9:۶۹ ال 


| M 


R 


duq ህ ህ MMM 





lure 0.50 Spectral Output with Two Passed Signals 


Figure 6.6 shows the overlap of both the curves of 


5.5 Note that the tvo signal spectral conponents 
is 


4 
Figure 
simply add and that their frequency relationship 


unchanged. 
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The effects of harnonically related signals vere 
demonstrated in the above data. These effects were used to 
aid in determining the frequency selectivity of the filter. 
This was done by usinq two input pulse trains which verə 
harmonically related: ۷.۱۷۷4۲ hertz and 8.88K hertz. The 
spectral components of the 8.88K hertz signal are 8.88K 
hertz, 17.76K hertz, 26.64K hertz, etc. and those of the 
4.44K hertz signal are 4.44K hertz, 8.88K hertz, 13.22K 
hertz, 17.76K hertz, etc. Thus, the line spectrum of the 
sum of these two signals will reveal that the spectral 
components of the 8.88K hertz signal should be larger than 
those of the 4.44K hertz ones because they are an addition 
of both the 8.88K hertz and 4.44K hertz frequency 
components. Figure 6.7a shows this effect on the linear 
spectrum plot. This spectrum is the addition of two signals 
(8.88K hertz and 4.44K hectz) and was used as the input to 
the filter. Figure 6.7b shows the filter output. In both 
pictures, the larger amplitude conponents “represent the 
summation of the 4.44K hertz and 8.88K hertz spectral 
components; and the smaller spectral lines between are only 


the 4.UUK hertz spectral lines. 
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Figure 6./a Input Spectral Components 
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Figure 6.7b Qutput Frequencv Spectrum with Two Input Signals 





The differences in spectral component amplitudes vas 
used in determining at vhat frequency the 8.88K hertz signal 
vas attenuated and not prəsent at the filter output. Then, 
by taking the difference of the passed signal frequency and 
the unpassed frequency,the filter selectivity could be 
determined. Figure 6.8 shows a comparison of the output of 
the 4.44K Hertz and 8.78K Hertz input signals. Note that 
the components of the 8.8K Hertz signal are deteriorating 
and decreasing. Decreasing the signal further to 8.59K 
Hertz as shown in Figure 6.9 shows that the original 8.88K 
Hertz is nearly gone with the exception of some 8.0K Hertz 
and 16K Hertz components. The filter output at 3.59K Hertz 
vas very distorted and the pulse train vas undetectable 
visibly on the oscilloscope. Thus, using this frequency as 
the limit of the filter"s passable frequency, results in a 
selectivity of 290 Hertz on either side of the passband 


center frequency. 
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Figure 6.9 Limit of Frequency Selectivity 
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C. PROBLEMS ENCOUNTERED IN DATA TAKING: 


One of the major problems encountered during data taking 
was the problem of the sampling frequencv clock drifting as 
much as 10K Hertz. Slight adjustment in the circuitry was 
continuously necessary. Also related to the clock problem 
was the problem of keeping the clocks in synchronization 
when the filter was connected in the second order design. 


However, once synchronised they generally remained so. 


The pulse generator was another source of difficulty in 
that th2 selected frequency for both the pulse Sorter and 
the SNR improvement demonstration would drift enough to 


cause the output to be distorted appreciably. 
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A. BACKGROUND 


The most common use of the recursive  canceller comb 
filter is as a cancellation circuit for a Moving Target 
Indicator (MTI) radar. In this type of system it is 
important to attenuate the unwanted signals (called clutter) 
from the radar return spectrum. This clutter results fron 
Stationary targets such as; mountains, buildings from the 
local surroundings, and moving targets such as the leaves on 
trees which are moving with the wind , the trees themselves 
swaying with the wind, ocean waves and swells, etc.. It is 
highly desirable that this narrow band of clutter signals be 
cancelled and the desired signal targets pass through the 
filter and eventually give a target presentation on the 
radar scope. The recursive canceller comb filter is ideal 
for this application. Since rt 15 periodic, a najoritv of 
the narrow frequency-banded clutter, which is periodic due 
to the pulse repition frequency (PRF) of the radar signals, 
will be cancelled. Also, the canceller frequency srectrun 
characteristic has a wide pass band and will pass the 
doppler shifted moving target. (unless it falls in the 
filters stop band which is then at the radars' blind speed). 
Compiete cancellation of the clutter will not OCCUR, 
however, since the stop band of the filter is generally not 
as wide as the clutter spectrum. If the clutter is a single 
spectral component, (stationary targets) then the filter 


Will cancel it. 
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B. CANCELLER FIITER DEMCNSTRATION 


The technique of cancelling an unwanted frequency and 
passing a wanted signal was demonstrated with the Feticon 
SAD-100 used previously in this study. The filter was 
connected as a first order comb filter. The clock frequency 
was reduced to its lovest possible value of 104.6 KHz. 
Although the filter is capable of operating at a clock 
frequency as low as z KHz, this frequency was not achievable 
without physically changing a capacitor end a variable 
resistor on the printed circuit board that the SAD-100 and 
its associated circuitry was mounted on. The Clock 
frequency of the filter was reduced to allow cancellation of 
as low a frequency as possible to simulate the previously 
discussed clutter Spectrum. However, at this clock 
frequency, the lowest frequency that could be canceilod vas 
DC and multiples of 1.09K Hertz. The frequency response of 


the filter at f =104.6 KHz is shown in figure 7.1. The stop 
C 
^A 


band center frequencies are at whole number multiples of 
1.09 KHz as seen in figure 7.1. For demonstration purposes 
the unwanted and desired signals were arbitrarily chosen to 
Ee 2 5582 and 13.1 Küz. Figure 2 depicts, in" the top 
oscilloscope trace, the addition oÉ the two signals at the 
filter input. “The bottom trace shovs the filter output of 
the 13.1 KHZ signal. Figure 7.3 shows the passage of 13.1 
KHz signal and a 1.5 KHz signal. This demonstrates that the 
filter will pass two simultanious signals and that if the 
unwanted signal is not at the stop band frequency, it will 


not be attenuated. 
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VIII: SUMMARY AND RECOMMENDATIONS 


Results of this thesis are summarized in three areas: 


Experiment, Theorv, and Applications. 
A, EXPERIMENTAL RESULTS 


Sampled analog recursive comb filters usinga first or 
second order CTD delay line have been studied experinentally 
in detail. The filter circuit used for this study vas the 
canonical circult universally used for the digital recursive 


filters. The feedforvard coefficients in this circuit Gar 


Bar and a) are restricted to a mutual relation of 1, +2, 1 


in order to achieve perfect separation (infinite 


attenuation) between neighboring comb. teeth. However, 
experimental results revealed that the nulls between comb 
teeth cannot be maintained. Deterioration of these nulls 
results with increasing frequency. To combat this problem, 
it was found that a sample and hold circuit was necessary at 
the input of the filter. The nulls of the filter frequency 
response improved in most cases giving nearly perfect 


attentuation between the neighboring comb teeth. 


The major portion of this thesis dealt with the 


integrator comb filter while a companion thesis (written by 


11] 
Arif  Ejas ) concentrated on the canceller comb features. 
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The Reticon-96 bit delay SAD-100 vas the device used to 
implement the comb filter. This device's characteristics 
were studied extensively and the measurements taken are 
contained in Chapter 4. This study revealed that sampled 
analog recursive comb filters deviated from the ideal comb 
characteristics in that: 

1. The peaks of the comb filter frequency response 
vary with frequency due to the Sin X/X characteristic of the 
Sample and hold circuit, and also with the frequency 
dependence of the numerator and denominator coefficients of 
the filter transfer function. 

ER A sample and hold circuit vas necessary at the 
front of the canonical filter circuit to improve the nulls 


in the frequency response. 
B. THEORETICAL ANALYSIS 


To model the sampled analog recursive comb filter from 
the existing digital recursive filter theory, it vas 


necessary to consider the folloving: 
1. The feedforvard coefficients Er a, and > must 


be related by the ratio of 1, +2, 1 respectively. This vas 
necessary since the coefficients will place a zero on the 
unit circle in the Z-plane at either Z = +1+j0 or Z = -1+30 , 


resulting in the desired nulls. 


zı 
2. The single delay Z used in digital recursive 


n 
filter theory, must be changed to Z where n is the number 


of delay stages in the device. This accounts for the n comb 


teeth from DC to the sampling frequency (f ). 
s 
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3. The filter coerficients are all dependent on 
frequency due to several causes: sampled and hold effect, 
frequency rolloff due tə limited bandwidth of electronics, 


and charge transfer inefficiency effect. 


In summary, the theoretical transfer function for a 


second order sampled analog recursive comb filter is: 


“ከ sen 
ao (f) +t a, (£)2 + ag (f)2 
HOZ RP رر‎ eeh 
1 + b (fS)Z " + b, (f) Z 
where 1 = 1, AR = +2, 85 = 1 at low frequency. 
L 


Using this modified theory, the theoretical analysis and 
experimental measurements are in close agreement. However, 
if proper frequency dependence of the filter coefficients 
can be obtained, the agreement between theory and 


experimental measurements will be improved. 


e. APPLICATIONS 


Three applications have been demonstrated, tvo for the 
integrator comb filters and one for the canceller type comb 


filters: 


1. Integrator comb filter 
a. Signal to noise ratio improvement 


b. Pulse sorting 


2. Canceller comb filter 
\ 


a. Deleting unwanted sinusoids 
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D. CONTINUATION OF THIS RESEARCH 


Suggestions for future vork are listed belov: 


1. 


Study the frequency dependence of the 
coefficients. This study vould involve more 
computer simulation and matching of actual 
data and the computer plots such as those 


contained in Chapter 3. 


Physically implementing the pulse sorter at the 
front end of a vide band receiver to evaluate 
its use under physical "real world" 


conditions. 


Implement the Reticon filter vith a chirped 
clock frequency and evaluate its performance. 
This type of device would have application 


in a chirped or jittered PRF radar. 


Determine the input-output phase relationship 
of the pulse sorter of Chapter 6 to determine 
if there is any phase distortion encountered 
by the filter. 


Determine the affects of the unwanted spectral 
components that fall in the filters pass 
bands on the desired pulse train spectral 
components for the pulse sorter discussed 


in Chapter 6. 


Theoretically compare actual data with the data 
of figure 5.26 based on area under the comb 


teeth versus the response attenuation area. 
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APPENDIX A 


PULSE TRAIN FOURIER ANALYSIS 


Pulse Width = t 


| | ۱ | 1 
=T 0 to T Sch ai 
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۱ TÈ -jX 
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Solving for the absolute value of oe 
lc. | = (Vt/T) ( SinX/X) 


Similarily solving for the DC value of the pulse train: 


T tort 
C =a = ۲ İsa = VT | var 
0 0 
ዐ to 
təl 
NN UTE 
0 to 
тәсис - 
0 ( ) V ( 0 0: 


Therefore: 


7— 50 
; (1/1) 
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APPENDIX B 


COMPUTER PROGRAM FOR PLOTTING HCZ) 


READ(5,200)A0,A1,A2,B1,B2 
200 FORMAT (7F6.3) 
PI=3. 141592654 
DEL=0. 1 
HMAX=0.0 
DO 5 ፲=1,961 
K=I-1 
IF (I-1) 2,2,4 
2 EGET: Ó 
GO TO 1 
4 GG=ABS (SIN (P1I*K*416.67/400000) /(PI*K*4U16.67/400000)) 
G CG 
1 CONTINUE 
C=COS (2*PI*K*DEL) 
S=SIN (2*PI*K*DEL) 
CC=COS (4*PI*K*DEL) 
SS=SIN (4*PI* K*DEL) 
XR=AO+A1*C +A 2*CC 
XI=A1*S+A2*SS 
YR=1+B 1*C+B2*CC 
YI=B1*S+B2*SS 
R= XR*YR+XI*YI 
RI=XR*YI-XI*YR 
0= ሂ፪ቾቾ 2+ፔ፲ *ቾ 2 
H(I)= (SORT (R**2+RI**2) /D) *GG 
IF (H(I) -HMAX) 6,3,0 
3 HM AX=H (I) 
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CONTINUE 
CONTINUE 

DO 9 ፲=1,961 
16 2 ۲-1 

F (I) =K*0.4 

A (I) =H (I) /HMAX 
CONTINUE 

CALL INIT 
CALL ХРЕМ(4) 
CALL YFRM (4) 
CALL NPTS (961) 
CALL PLOT (F, A) 
CALL PAUSE 
CALL FIN 

EN D 
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APPENDIX C 


EXPERIMENTAL DATA FOR SNR GRAPHS 


Coeff. Noise Volts Signal Volts SNR 
(RMS) (RMS) (DB) 
b ۷ V V ሃ in Out Improv. 
1 0 i 0 1 

+0.8 0.061 0.038 0.036 0.021 -5.36 -4.58 40.78 
+0.7 0.059 0.405 0. 5 0.024 -4.51 -4.55 -0.04 
+0.6 0.056 02345 02035 0.202 -5.37 -4.00 +1.37 
+0.5 0.056 00059 0.036 0.226 “5757155521: ٦2 
40.2 0.055 0.037 0.039 0.023 23.90 -2.99 3 

0:0 0.063 07039 0.039 0.022 -5.18 -4.17 +1.02 
-0.2 0.068 0.037 0.043 0.020 -5.37 -3.98 41.39 
-055 0.079 0 19 0,062 07022 SSMO -2.13 +F2795 
206 0.097 0.038 0.083 0.021 -5.25 1.31 +3794 
“057 0.118 0.042 0. 114 0.023 =5.21 -0.30 +4.91 
-0.8 0.114 0.038 0.195 0.021 -5.36 +2.63 +7.99 


*Data plotted as dots with squares around them 


ia Figure 5.26. 
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Coeff. 


40.8 
40.7 
10.6 
40.5 
+0.4 
2 
572 
+0.1 

0.0 
-0.1 
UM 
053 
-0.4 
о 
x: 
057 
-0.8 
-0.9 


Noise Volts 


(RMS) 
V V 
0 E 

0.126 0.082 
0.115 0.082 
0.119 0.085 
0.120 0.085 
0.123 0. 085 
0.125 0.085 
0.125 0.085 
0.135 0.085 
0.134 0. 085 
0.135 0. 085 
0.150 0. 085 
0.152 0. 085 
0.162 0.085 
0.196 0. 085 
0.219 0. 085 
0.244 0. 082 
0.320 0.082 
0.575 0.085 


Signal Volts 


(RMS) 
۲ Vo 
0 pk 

0.066 0.044 
0.060 0.044 
0.060 0.044 
0.063 0.045 
0.070 0.044 
0.077 0. 055 
0.076 0.045 
0.088 09045 
0.094 0.045 
0.088 0.045 
0.099 0.045 
0.098 0.045 
0.125 0.044 
0.181 0.045 
0.218 0.044 
0.295 0.044 
0.425 0.044 
1.300 0.044 


In 


2-39 
-5.43 
məz 
-5498 
EA 
-5 060 
=ጋ و‎ 
-00 
د وت‎ 
- 2 
> اه‎ 
-5.60 
xo50 
S58 
-5 ۶2 
-5 3 
2-5393 
S2 


SNR 
(DB) 


Out 


=5. 68 
- 97 
-52 6 
obl 
-4.91 
74.26 
=i. 38 
207 
2.02 
290702 
1 
-3.81 
22925 
-0.69 
-0.04 
41.65 
+2.46 
17:09 


*Data plotted as dots with triangles around then 


11 5210503126 5.6. 
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12020۷ 


-0.15 
-0.24 
=0 24 
-0.03 
EL 
+1. 
+0 
41.84 
12.54 
+1.90 
42.91 
+1.79 
+3.45 
+4.89 
555: 43 
47.08 
5 
412.28 








Coeff. Noise Volts Signal Volts SNR 

(RMS) (RHS) (ƏB) 
b V v V V In Out  Inprov. 

1 0 ı 0 1 

+0,7 0.204 0. 150 0.038 0.023 -16.5 -14.6 56 
+0.6 0.206 0.150 0.037 0.023 -16.5 -14.9 +1.56 
+0.5 05207 0.150 0.037 0.023 -16.5 -14.9 +1.59 
+0. 4 0.212 0.450 0.039 0.023 -16.5 -14.6 41.93 
+0.3 0.220 ) 0 0. 2 0.023 -16.5 ~14.4 +2109 
5 0.228 0. 150 0.043 0.023 =16.5 -14.5 +1.99 
40.1 0.236 0.150 0.7 0.023 -150:720110 0 ус 
0.0 0.235 0. 0 0. 5 07023 -16.5 -14.3 +2.18 
=0.1 Ur 2 45 0. 150 0.047 0.023 -16.5 -14.0 +2.46 
02 0.242 0.150 0.049 0.023 -16.5 -13.9 +2.61 
-0.3 0.261 0.150 0.056 0.023 219 uzun + m 
US 0. 0 0. 0 0.063 0.023 216.5 1229 69 
S 02312 0.150 0.073 0.023 =16.5 =-12.7 +3.80 
-06 07353 0.150 0.094 0.023 2-160959 -11:5 +4.98 
UE 0.410 0. 150 0.123 0.023 -16.5 -19.5 46.02 
IUE 0.525 0.150 05197 0.023 =16.5 051670 
E 0.905 0. 150 0.480 0.023 >|ርር ያ ፡ምሜን َ٘٘٘,[ 9 


*Data plotted as dots with circles around them 


in Figure 5.26. 
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